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Abstract 

We present an analytical method of calculating the mean first-passage times (MFPTs) for the 
magnetic moment of a uniaxial nanoparticle which is driven by a rapidly rotating, circularly po- 
larized magnetic field and interacts with a heat bath. The method is based on the solution of the 
equation for the MFPT derived from the two-dimensional backward Fokker-Planck equation in the 
rotating frame. We solve these equations in the high-frequency limit and perform precise, numer- 
ical simulations which verify the analytical findings. The results are used for the description of 
the rates of escape from the metastable domains which in turn determine the magnetic relaxation 
dynamics. A main finding is that the presence of a rotating field can cause a drastic decrease of 
the relaxation time and a strong magnetization of the nanoparticle system. The resulting station- 
ary magnetization along the direction of the easy axis is compared with the mean magnetization 
following from the stationary solution of the Fokker-Planck equation. 

PACS numbers: 75.60.Jk, 76.60.Es, 75.50.Tt, 05.40.-a 
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I. INTRODUCTION 



The problem of finding tlie statistical characteristics of the first, the biggest, and the like, 
for sample paths of a stochastic process frequently occurs in physics, biology, economics, and 
other sciences .-i2i>^ A particularly prominent identifier is the mean first-passage time (MFPT), 
i.e., the average value of the random times taken by a random walker that starts out from 
some initial state to reach another prescribed state in state space for the first time. This 
quantity describes a large variety of noise induced effects such as activation rates or reaction 
rates, the lifetime of metastable states, the extinction of populations, or the extreme events 
in financial time series, to name only a few. Unfortunately, the class of stochastic processes 
for which the MFPT can be calculated explicitly is rather restricted. In fact, the most general 
analytical results were obtained for one-dimensional discrete or continuous Markov processes 
that are homogeneous in time.-*^'^'^ However, these conditions of one- dimensionality and 
time-homogeneity often represent an oversimplification. In particular, Markov processes 
that describe time- dependent systems are usually not homogeneous. Prominent examples 
are Brownian motors and ratchet-like stochastic systems,-»2, as well as systems exhibiting 
Stochastic Resonance.— Although a few advanced, approximate methods for the analysis 
of periodically driven Markovian systems are available^'^^'i^ii^'i^ii^'i^ii^ the development of 
new approaches for calculating the MFPTs in such systems still presents an important 
challenge. 

In this paper, we develop an analytical approach to the two-dimensional MFPT problem 
for a magnetic moment of a ferromagnetic nanoparticle driven by a magnetic field which 
rapidly rotates in the plane perpendicular to the easy axis of magnetization (up-down axis). 
The natural precession of the magnetic moment always occurs in the counterclockwise di- 
rection (when viewed from above). Therefore, its deterministic dynamics in the up and 
down states differ.— For this reason the stochastic dynamics and thus the MFPTs in these 
states are different as well. In turn, the difference in the MFPTs can drastically change 
the magnetic properties of systems composed of nanoparticles. This was explicitly shown 
in the resonant case, when the driving frequency coincides with the Larmor frequency of 
precession.— With this work we present a detailed analysis of this phenomenon in the case 



of a rapidly rotating magnetic field, which was briefly presented in Ref. 



22|, and apply the 



MFPTs for describing the thermally activated magnetic relaxation in such systems. 
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The paper is structured as follows. In Sec. II, we present the model and introduce the 
stochastic Landau-Lifshitz equation together with the corresponding forward and backward 
Fokker-Planck equations in the rotating frame. In Sec. Ill, we derive the general two- 
dimensional equations that define the MFPTs for the driven magnetic moment of a uniaxial 
nanoparticle in the up and down states. The analytical solution of these equations in the 
case of a rapidly rotating magnetic field is carried out in Sec. IV. In the same section, to 
verify our method, we solve the effective stochastic Landau-Lifshitz equations and calculate 
the MFPTs numerically. Some apphcations of the obtained results are presented in Sec. V. 
Here we study the features of magnetic relaxation and steady-state magnetization induced 
by a rapidly rotating field in nanoparticle systems. We summarize and discuss our findings 
in Sec. VI. 



II. MODEL AND BASIC EQUATIONS IN THE ROTATING FRAME 

Let us consider a simple model of a uniaxial ferromagnetic nanoparticle within which 
the nanoparticle is characterized only by the anisotropy field and the magnetic moment 
m{t) of fixed length |m(t)| = m. This model is relevant for nanoparticles whose sizes do not 
exceed the exchange length, i.e., the length scale below which the exchange interaction is 
predominant. We assume also that perpendicular to the easy axis of magnetization, which 
we choose as the z axis of a Cartesian coordinate system xyz, a circularly polarized magnetic 
field h{t) — h{cos Lot, p sin Lot, 0) is applied. Here, h — \h{t)\ is the field amplitude, lo is the 
angular field frequency, and p = — 1 or for clockwise and counter-clockwise rotation of 
h(t), respectively. 

We take into account the influence of a heat bath by means of damping and the presence 
of a thermal, Gaussian distributed magnetic field n(t), possessing zero mean and the white 
noise correlations {na{ti)ni3{t2)) — 2ASai3S{t2 — ti). Here, na{t) (a,/? = x,y,z) are the 
Cartesian components of n(t) , A is the intensity of the thermal field, is the Kronecker 
symbol, S{t) is the Dirac 5 function, and the angular brackets denote averaging with respect 
to the sample paths of n(t). In this case the dynamics of m(i) is Markovian and can be 
described by the stochastic Landau-Lifshitz equation 

— m = -7m X (Heff + n) m X (m X Heff), (2.1) 
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where 7 > denotes the gyromagnetic ratio, A(> 0) is the dimensionless damping parame- 
ter, Heff = —dW/dm is the effective magnetic field acting on m(t), W denotes the magnetic 
energy of a nanoparticle, and the cross (x) indicates the vector product. The stochastic 
Landau-Lifshitz equation in the form of fl2.ll) was introduced by Kubo and Hashitsume.— 
Another form of this equation (with Gilbert's relaxation term^^) was employed by Brown 
in his well-known paper.— We note, however, that although the solutions of these equations 
for a given realization of n{t) are generally different, their statistical characteristics are the 
same^ (to within a renormalization factor for A). At present, both forms of the stochastic 
Landau-Lifshitz equation are used equally often. Although in the Langevin equation (12. ip 
the noise enters in a multiplicative manner, the resulting process actually is independent of 
the stochastic calculus if restricted to the sphere m^(t) = m^. This is so, because the corre- 
sponding, noise induced drift terms are perpendicular to the sphere for any interpretation of 
the stochastic differential equation (I2.ip .— '^'^'^ For the sake of definiteness, upon perform- 
ing nonlinear transformations of this stochastic differential equation we shall consistently 
employ the Stratonovich interpretation.— We next perform such a nonlinear transformation 
by considering the dynamics of m{t) in terms of spherical coordinates in a rotating frame. 

Specifically, using the respective polar and azimuthal angles 6 and ip of the magnetic 
moment, m(t) = m(sin ^ cosy), sin ^ sin y), cos ^), the magnetic energy W emerges as 

= - mHa sin^ 9 — mh sin 9 cos{{p — put). (2.2) 

The energy W depends on (f and t only through the single variable ip = (f — put. Therefore, 
it is advantageous to introduce a rotating Cartesian coordinate system x'y'z' (see also a 
similar description of a noisy, periodically driven Van der Pol oscillator in Ref. 3]), in 
which h(t) = h{l, 0, 0) and the azimuthal angle ip is changed by ip. According to Eq. (12.11) . 
in this coordinate system the equations for 9 and ip can be written in the dimensionless form 
as follows: 

9 = u{9,iP) + U0,i^,r), 
iP = v{9,iP)- pn + ^49,iP,T). (2.3) 

Here, an overdot denotes the derivative with respect to the dimensionless time r = Urt 
with the Larmor frequency cu^ = jHa, and Q = uj/ujr is the corresponding dimensionless 
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frequency of the driving field. The functions u{9,'4)) and v{6,iIj) result as 

«(»,</■) 



sm^ 9 \ do dip 



(2.4) 



with W = W{0,ip) /mHa denoting the dimensionless magnetic energy. The stochastic forces 
ie,i){0,4',T) are determined as 



= —n^i (r) sin ip + Uyi (r) cos ip, 
= (r) — [hx' (r) cos + ^y' (t) sin ip] cot 6*, 



(2.5) 



where nx>{T) = nx{T) cos Qt + pfiylr) sin Qt, fiyiij) = fiylr) cos Qt ~ phx{T) sin Qr, and 
nai^) = ^^{t / ujr) / Ha are the components of the reduced thermal magnetic field. Using 
the statistical characteristics of na{t), for these components we readily obtain {hair)) = 
and (nQ,(ri)n^(r2)) = 2A6a/36{T2 — ti). Here, A = Aj/Ha is the dimensionless intensity of 
the reduced thermal field. With the help of the Sutherland- Einstein relation,— i.e., A = 
XksT/'-fm wherein fc^ denotes the Boltzmann constant and T is the absolute temperature, 
it can be written also in the form A = X/2a, where 



a = mHa/2kBT 



(2.6) 



is the anisotropy barrier height in the units of the thermal energy ksT. We note that in the 
purely deterministic case, when A = 0, some important features of the solution of Eqs. (12. 3p 
were studied in the context of the nonlinear dynamics of m(t) and its stability. ^-i^ 

Next we introduce the conditional probability density P = P{6,tp,T\6' ,tp' ,t') (r > r') 
which presents the most important statistical characteristic of the solution of Eqs. (12. 3p . Us- 
ing the well-known connection between a set of stochastic differential equations and the corre- 



sponding Fokker-Planck equation, see, e.g., Refs. 
forward Fokker-Planck equation 



29 



37l | , we obtain the two-dimensional 



2a 



d^P 1 

+ 



d'^p 



(9^2 sin'edip'^ 



d_ 
86 



— cote + u(9,i!) 

2a 



P 



(2.7) 
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and the corresponding backward Fokker-Planck equation 



A \d^P 1 ^^P^ 



dP 



2a [d9'^ sin^e'dtp'^ 




dP 



(2.8) 



Notably, this two-dimensional Fokker-Planck dynamics does not obey a detailed balance 
symmetry,— >^ if the driving frequency Q does not vanish. In the model under consideration 
we have W = (l/2)sin^^ — /i sin cos'?/', where h = h/Ha- Therefore 



The probability density P must satisfy the equal-time condition P\t=t' = S{9—9')6{ip~4'') 
and appropriate boundary conditions, as implied by the physical context. Moreover, in spite 
of the singularities in Eqs. (12. 7p and (12. 8p at 6, 6' = 0,71 (which are a consequence of the use 
of the spherical coordinate system), the probability density P must be a regular function 
also at these points. In addition, if not excluded by the boundary conditions, P must be 
properly normalized, i.e., Jq^ dip d6 P = 1. We note also that the forward and backward 
Fokker-Planck equations (12.71) and (12.81) are equivalent; the difference between them is which 
set of variables, 0,ip,T or 6',ip',T', is held fixed. Due to this difference, the former is more 
convenient for studying the statistical properties of the magnetic moment m(t) as functions 
of the evolving time t, while the latter one is more appropriate in studying the first-passage 
time statistics for m(t). 

Based on the Fokker-Planck equation (12.71) one can determine stochastically equivalent 
Langevin equations, reading 



where rjei'^) and i]^{t) denote two independent Gaussian white noise sources with zero 
mean and white noise correlations {f]i{'^)Vji'^')) = ^ij^i''' " t')? wherein i,j = 9,ip. In spite 
of the multiplicative nature involving the noise ?7^(t) in the second equation, the resulting 
stochastic dynamics possesses a vanishing noise-induced drift and thus is again independent 



= —A sin 6* cos 9 + h{X cos 9 cos ip — sin ip) 



^ ~ cos 6' cos -j/^ -|- A sm 

[9, %p) = cos 9 — h . 

sm 9 



(2.9) 




(2.10) 
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of the employed stochastic calculus.—*^ As a basis for numerical investigations the latter 
Langevin equations provide a more convenient starting point than Eqs. (12. 3p : this is so 
because they require the simulation of only two, rather than three independent Gaussian 
white noises. 

III. GENERAL EQUATIONS FOR THE MFPTs 

In the absence of the random magnetic field, i.e., for n(t) = 0, or equivalently A = 0, the 
motion of the magnetic moment follows the deterministic Landau-Lifshitz equations. In the 
rotating frame the resulting deterministic dynamical system is not explicitly dependent on 
time and given in terms of two degrees of freedom. For A > this constitutes a dissipative 
system, which can only perform regular motion approaching fixed points or limit cycles in 
the asymptotic limit of large times. In the present paper we are mainly interested in values 
of the parameters A, h and Vt for which the motion is bistable, i.e., the asymptotic motion 
does lead to either of two attractors depending on the initial condition. These attractors 
are denoted as up and down states and labelled by a = +1 and cr = — 1, respectively. The 
dynamics generates a partition of the state space into two domains of attraction, containing 
either the up or the down state. The common boundary between the domains of attraction 
is formed by the separatrix. 

In the presence of a random magnetic field the separatrix is no longer an impenetrable 
border, and transitions between the two domains of attraction may occur. For small random 
fields, corresponding to large values of a, these transitions are rare and can be characterized 
by transition rates. It is now tempting to determine the rate between a state a and the 
opposite state by the MFPT to the separatrix, T^sep, i-e., by the statistical average of the 
stochastic first-passage times of trajectories that start at the attractor a and reach the 
separatrix for the first time. In the asymptotic limit of vanishing noise a trajectory visiting 
the separatrix will go to either side with equal probability and the rate is given by the inverse 
of twice the mean first-passage time. For finite noise the transition from the separatrix to 
the two sides may differ from each other,— whereby the precise value of the resulting bias 
in general is difficult to quantify. In order to be more flexible, we consider the mean first- 
passage times to the boundaries of two regions = {9,ip\0 < < (f)^i{ip),0 < ip < 27r}, 
and R-i = {6,tlj\(f)-i{ilj) < 6 < tt,0 < ip < 27r}, each containing one attractor. The 
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boundaries 4>a{'ip) can be chosen as the separatrix or as any other curve between the two 
attractors. For a convenient choice of the boundaries 0cr(V') 'we refer the reader to the next 
section. We note that the regions are stationary in the rotating frame but may move in 
the rest frame. 

In order to determine the first-passage times of the regions R„, re-crossings of the re- 
spective boundaries dR^ = {0,ip\9 = (j)a{i'),0 < ip < 27r} must be suppressed.- This is 
conveniently achieved by imposing absorbing boundary conditions on the conditional proba- 
bility density PaiO-, ipiT — t'\9' , ip') obeying the backward Fokker-Planck equation (12.81) . i.e., 
we require Pa{9,ip,T — T'\9',tp') = for {6',ip') G dR^. Here, we used that in the rotating 
frame the process in the stationary regions R^j is time-homogeneous with 

p„{e, i),r- T'\e', i)') = p^{e, r - r'\e', o) 

= p^{e,ij,T\e',ij',T'). (3.1) 

If the magnetic moment starts out at the time r' at the position {6', G Ra it will unin- 
terruptedly stays within the initial region R„ until a time r with a probability Q{6', ip'; t—t'). 
This probability can be expressed as the integral of the conditional probability density over 
all states in R^, i.e., 

Q^{e\^'-T-T')= [ d^dep„{e,^,T-T'\e'^'). (3.2) 



Ra 

The probability Qcr{0',il''',T — t') is a solution of the backward equation with the absorb- 
ing boundary conditions Qa{0',ip'',T — t') = for {6',iIj') G dR^ and the initial condi- 
tion Qcr{0' = 1- Integrating Qcr{6' ,iIj';t — r') over all positive (dimensionless) times 
u = r — t', one obtains an expression for the (dimensionless) MFPT of the form 

POO 

%{e\^') = / duQ„{e\^'-u). (3.3) 

JO 



This MFPT is the solution of the backward equation 



2^ 



dO''^ sin^ 9' dip''^ 



+ 



— cot 9' + u(9', f 
2a 



89' 



+ K^',^')-Pfi]^ = -1 (3.4) 
with the absorbing boundary conditions 

%{9\ = for {9', ij') G dR^. (3.5) 



Eq. (13.41) was derived in Ref. [38|] for the undriven case with ri = 0. 

Because u{9',ip'), v{9',ip') and also Xj{9',iIj') are periodic functions of it is conve- 
nient to decompose these functions into their average and periodically varying parts in ip': 

Here, ui{9',ip') = vi{9',ip') = Sa{0',ip') = 0, the overbar denotes an average over ip', i.e., 
00 = (l/27r) /o^'^ di;'{-) and, according to (K9\f . 

u{e') = -A sin 9' cos 6', v{e') = cos 6', 



Ui [6', i/j') = h{X cos 6' cos ip' — sin ip') , 

~ cos 6^' cos ^z^' + A sin ?/^' 

^ , ¥) = -h . 

smt'' 



(3.6) 



Using these decompositions, we find from Eq. (13.41) coupled equations, see also in Ref. |22| . 

for the average part 



2a 



r + cot e' 



d9' 



dip' 



(3.7) 



and for the periodic part 



2a 



09'^ sin^ 9' di)'^ 
+ (w - pfi) 



1 d^S dS^ 

+ cot 9 



89' 



dTn 



ur 



d9' 



u 



09' 



dScr dSa- dSa- 



Ui- 



di,' -09' = ^'-'^ 

We emphasize that these equations are fully exact, i.e., they follow from the stationary 
backward Fokker-Planck equation being in the rotating frame. 



IV. RAPIDLY ROTATING FIELD 



A. Analytical analysis for the MFPT 



In the case of a rapidly rotating magnetic field, i.e., when the condition h/Q <^ 1 holds 
(note that h need not be small), Eq. (13. 8p can be essentially simplified. The reason is that 
the function and its derivatives tend to zero as h/Q — > 0. Using these conditions and 
taking into account that for large frequencies the term QdS^/dip' is of the order Q^, we 
obtain from Eq. (13. Sp the approximation: 



OS, 
dip' 



dTa_ 
d9' 



0. 



(4.1) 



According to (13.61) . the solution of Eq. (14.11) . which satisfies the condition = 0, is given 

by 

5, = p-(Acos^'sin^' + cos^')^- (4.2) 

This solution self-consistently conforms with the assumptions made above. Using (14. 2 p and 
dMD we find _ 

where 

ks = -ph^n. (4.4) 
With these results, Eq. (13. 7p reduces to the form 

A (^^T A ~ rIT' 

According to this equation, a magnetic field rapidly rotating in the plane perpendicular 
to the easy axis of the nanoparticle acts on the nanoparticle's magnetic moment precisely 
as a static effective magnetic field hes (in units of the anisotropy field Ha) which is applied 
along the easy axis. The direction of this field and the direction of the field rotation follow 
the left-hand rule and its value can be large enough to produce observable effects. Next we 
assume that \hes\ < 1; otherwise only one state, cr = +lor(j = — 1, is stable. 

It is not difficult to show that the general solution of Eq. (14. 5p contains both regular and 
singular parts. The singular part arises solely from the use of a spherical coordinate system 
and has no physical meaning. It exhibits a logarithmic singularity at 6' = 7r(l — a)/2 and, 
as a consequence, its derivative diverges. On the contrary, the regular part has a vanishing 
derivative at this point. Therefore, in order to exclude the contribution of the singular part, 
the solution of Eq. (14. 5 p must satisfy the regularity condition {dTa-/d9')\oi=T,(^i^„)/2 = 0. 
We note in this context that the regularity condition corresponds to the situation when a 
reflecting barrier is placed at the point 6' = 7r(l — o")/2. 

In order that the high frequency approximation for the MFPTs can consistently be per- 
formed, the ip dependence of the boundary curves (paii^) niust be chosen conveniently. 
Assuming that (paii^') = 0o- + (piaii^') with (?/;') = and cpi^^ip') ~ h/Q, the absorb- 
ing boundary condition (13. 5p in the linear approximation in h/Q leads to the relations 
7a((pa) = and (pia{i^'){dT^/d9')\g,^-^^+Sa(^„,ip') = 0. Using (H^]) and the latter relation. 
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we find for tlie absorbing boundary tlie explicit form 

(paiip) = <Pa- p^(Acos0^sinz/' + cosip). (4.6) 

Next, solving Eq. (14.51) with the specified regularity and boundary conditions, 
{dT^/d6')\e'=n{i-a)/2 = and T^(0^) = 0, we obtain 

TAO') = ^ / _ dx—^ — / dy e^^y+'^^^^\ (4.7) 

A J cos (f>^ i — X Jx 

where 9' G [0,0+^] if cr = +1, and 9' G [0_]^,7r] if cr = —1. The angles 0^ can be chosen 
depending on physical situation. For high potential barriers, i.e., a ^ 1, the magnetic 
moment predominantly dwells in the vicinity of either of two equilibrium states at 6* = 
and 9 = 71. In this case the transition times between the states a and —a exceed by far the 
relaxation times towards these states. Therefore, in dimensional units the averaged MFPT 
Tfj{9') = T(j{9')/ur representing the transition time from one state a to the opposite state 
—a only weakly depends on the precise location of the initial magnetization, as long as 9' lies 
within the domain of attraction of the considered state a. Accordingly, the precise location 
of the absorbing boundary (p^ practically has no effect on T„{9') if it is located well beyond 
the separatrix which divides the state space into domains of attraction of the up and down 
magnetization. Under these conditions, we find from (14. 7p in leading order in a: 

^ _ 7^ _ 1 exp [a(l + aKsf] ^^^^ 

As it follows upon inspection from (14. 2 p and (14.81) . in the high-frequency limit the periodic 
part of T^{9', ip') can be neglected, i.e., T^{9', ijj') ^ T^. 



If \hes\ ^ 1 then (14. 8 p yields T^- = Tq exp{a2ahes) , where To = {1 / XuJr) ^ it / a expa is the 
MFPT at h = 0. According to this formula, a rapidly rotating magnetic field increases the 
MFPT for the magnetic moment in the state a = —p {ahcs > 0) and lowers this MFPT 
for the magnetic moment in the state a = +p (cr/iefr < 0). This difference in the MFPTs 
arises from the natural precession of the nanoparticle magnetic moments, which occurs in 
the counter-clockwise direction, if viewed from above. As a consequence, the statistical 
behavior of the up and down magnetic moments in the magnetic field rotating in a fixed 
direction is different. 

Another choice of the position for the absorbing boundary can be made right on the 
separatrix itself. For the averaged one dimensional dynamics of the azimuthal angle it 
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corresponds to the 6 value where the deterministic part of the drift of the reduced backward 
equation (14.51) assumes an unstable fixed point, i.e., for cos6' = —hes, or, equivalently, to the 
maximum of the effective magnetic energy of the nanoparticle, Wefr(^) = ''^-f^a[(l/2) sin^ — 
hcs cos 6]. Accordingly, the averaged dimensionless MFPT from an initial angle 6' to the 
separatrix reads 

T^so^e') = ^ I dx—^ — / rfye"(^+'^-)^ (4.9) 

A J —h^ff i — X Jx 

If a ^ 1 and cos 6*' is not too close to —h^s then To-scp(^') only weakly depends on 6' and 
Tcrscp Tcr/2 as a ^ oo. Thus, the ratio T^scp/Ta between the MFPT to the separatrix 
and to an angle (f)^ which is well beyond the separatrix converges to the value 1/2 if a oo] 
however, the larger the rotating field amplitude h the slower is the convergence, cf. in Fig. 1. 
We note also that for large but finite a the conditions TaseplTa > 1/2 and T^sep/Ta < 1/2 
hold for ap = +1 and ap = —1, respectively. The reason is that the effective magnetic energy 
Wcs{9) has different slopes from the left and from the right of the separatrix disposed at 
9 = arccos(— /leff)- 

B. Numerical simulations 

In order to examine the analytical results developed for calculating the MFPTs for a 
magnetic moment which is driven by a rapidly varying circularly polarized magnetic field, 
we performed numerical simulations of the full two-dimensional Langevin equations (12.101) . 
For any given set of parameter values p, cr, Q, a, A, h four groups of 10^ trajectories were run 
using a stochastic vector Euler algorithm.-i^ All trajectories of a simulation were initialized 
with the same values for 6 and ip. In the vicinity of the coordinate singularity at 6 = 
a reflecting boundary was located at = 0.027r, together with an absorbing boundary at 
6 = O.Stt, which is located well beyond the separatrix of the corresponding deterministic 
system. The step width was chosen such that the increments in 6 were less than vr/lOO and 
the increments in ip less than vr/10, respectively. The precise computation of the MFPTs 
of this system using the Langevin equations (12.101) requires to compute the arrival of all 
trajectories at the absorbing boundary. In practice this method is unfeasible for all but 
the smallest values of the anisotropy barrier height a; this is so because a given trajectory 
can take much longer than the MFPT to arrive at the absorbing boundary. These events 
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are rare, but contribute significantly to the MFPT, and hence a sufficiently large number 
of these events needs to be simulated to arrive at a reliable statistics. However, assuming 
that the ffist-passage time distribution is exponentially distributed with a rate parameter 
l/Tcr it is possible to determine the MFPT approximately by fitting the tail of the simulated 
first-passage time distribution to an exponential. Although this assumption can be justified 
for large barrier heights a ^ 1 since then the relaxation time of the system is much shorter 
than the MFPT, it clearly is expected to fail for a ^ 1. Note, that a enters the exponent of 
the expression for the MFPT equation (14.71) and therefore, the large barrier limit is already 
obtained for moderately large values of a. Using this assumption allows to simulate the 
Langevin equations fl2.10p for each trajectory up to a fixed maximal time at which on average 
a considerable fraction of all trajectories, but not all, have crossed the absorbing boundary. 
To be definite, we took this time to be two thirds of the theoretical mean first-passage time 
calculated from equation (14.71) . With this choice, roughly half of all trajectories arrived 
at the absorbing boundary. The number of absorbed transitions was stored as a function 
of time. From each of the resulting four data sets the rate 1/Ta- was determined by an 
exponential fit and thereby the MFPT was estimated. From these four values an average 
value and a standard deviation was determined. The so obtained numerical findings compare 
most favorably with our theoretical predictions, cf. Figs. 2 and 3. 

Fig. 2 depicts the dependence of the MFPT as a function of the anisotropy barrier height 
a for two different driving angular frequencies i7 = 5 and = 10. The values of the other 
parameters are p = +1, a = +1, h = 1, A = 0.1. Depicted is the natural logarithm of the 
dimensionless MFPT T+i. It can clearly be seen that the high-frequency predictions from 
the equations (14. 7p and (14.81) approach the results of the numerical simulation from above 
as Q is increased. This can be understood intuitively because the effective dynamic barrier 
for escape is increased as Q is increased. As expected, the agreement decreases in quality 
upon lowering the angular driving frequency; this fact is corroborated with the numerical 
results for Q = 5, cf. the thick dashed line versus the numerical data points. 

Fig. 3 shows the dependence of the MFPT as a function of h, A and a at constant a = 5 
and p = +1. As predicted by equations (14.71) and (14.81) . T„ can be seen to be proportional to 
1/A. The up and the down states {a = +1 and a = —1) are clearly inequivalent for non-zero 
h. This is in perfect agreement with the prediction of the theory, as can be deduced from the 
approximate equation (14. 8 p in which a = +1 (cr = — 1) decreases (increases) the exponent. 
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since hes < for p = +1. 



V. RELAXATION OF THE MAGNETIZATION 
A. Relaxation law at high anisotropy barrier 

The MFPTs provide important characteristics for the magnetic dynamics of nanopar- 
ticles. If their states a are used for information storage then the average times during which 
the information in these states is kept safely must be considerably shorter than To-. The de- 
pendence of Tcr on the characteristics of the rotating field gives a possibility to intentionally 
change the relative stability of the up and down states. 

The transition times T„ from the state a to the state —a also determine the thermally 
activated magnetic relaxation in a system composed of uniaxial nanoparticles whose easy 
axes are perpendicular to the plane of field rotation. If a ^ 1 and the precession angle^^ 



9, 



' + h (5.1) 



^ (1 - apny + A2 

of the magnetic moment in the state a is small, i.e., 6l 1, the reduced magnetization of 
this system can be defined as fi{t) = [A^+i(t) — N^i{t)]/N, where A^o-(t) denotes the number 
of nanoparticles in the state a at time t, and the whole number of particles. Taking into 
account that N+i{t) + N^i{t) = N, this definition yields fi{t) = 2N+i{t)/N. Next, since the 
rate is the probability of reorientation of the magnetic moment from the state a to the 
state —a per unit time, we have N^it) = N_^{t) /T — N^{t)/T^ and thus the equation for 
/i(t) takes the form 

+ + (5.2) 

Its solution with the initial condition /i(0) = 1 is given by 

/i(t) = (1 - Poo) exp(-t/trel) + yUoo, (5.3) 

where trei = T+iT^i/ {T^i + T_i) is the relaxation time and 

Poo = ^"^^ , (5.4) 
J +1 + i -1 

is the steady-state magnetization of the nanoparticle system induced by a rapidly rotating 
field. 
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A simple analysis shows that a rapidly rotating magnetic field causes a decrease of the 
relaxation time and magnetizes the nanoparticle system along the easy axis of magnetization 
(trei < ^rei and /ioo 7^ if /i 7^ 0). This conclusion is not trivial because the rotating field has 
no component in the direction of the induced magnetization. The value of this magnetization 
grows with decreasing temperature and its sign is determined by the direction of the magnetic 
field rotation, i.e., sgn/ioo = —p- In particular, if \hcs\ ^ 1 then t^ei = {To/2) cosh~^(2a/iefj) 
and yUoo = tanh(2a/ieff)- These results indicate that even a weak rotating field can drastically 
decrease the relaxation time and strongly magnetize the nanoparticle system if the temper- 
ature is low enough that a|/ieff| ^ 1- We note also that Eq. (15.41) overestimates the absolute 
value of the magnetization because the magnetic moment is less than m in the up state and 
larger than —m in the down state. 



B. Steady-state magnetization at high frequencies 



In the case of a rapidly rotating field we are able to calculate the steady-state magne- 
tization /ioo for an arbitrary anisotropy barrier. To this end, we introduce the stationary 
probability density Pgt = PstiP,'^) which, according to Eq. (12. 7p . satisfies the stationary (in 
the rotating frame) Fokker-Planck equation 



2a 



d9^ 
d 



+ 



d_ 



2a 



cot 9 + u{9,tp) 



(5.5) 



being properly normalized, Jq^ dip Jq dOPst = 1- Using the decompositions Pst = Pst(^) + 
Pi{9, i:) (Pi = 0), u{9, ij) = u{e) + ui{e, ij) (Hi = 0), and v{d, ^j) = v{9) + Vi{9, ^j) (yi = 0), 
where the overbar denotes averaging over ijj, i.e., (■) = (l/27r) Jq^ dilj{-), we obtain from 
Eq. (15.51) coupled equations for Pst 



St 



2a dO^ 



d_ 

de 



2a 



cot 6 + u 



- d 

Pst - uiPi 

ad 







(5.6) 



and Pi 



A ra^Pi 1 d^Pi 

+ 



2a 



09^ sin' 6 dip^ 



d_ 

de 



— cote + u(e,i!) 

2a 



Pi 



9 ^ d 

^«tPst + ^«lPl 
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(5.7) 



Assuming that Pi ~ f2 ^ as i7 ^ oo, Eq. (15. 7p reduces in the high-frequency hmit to 

(5.8) 



Since Vi = — /i(cos 6 cos ip + X sin ip) / sin 6 and Ui = h{X cos ^ cos ip — sin ■?/') , the last equation 
has the solution 

h, . Pst 



Evaluating the average 



Pi = — p— (cos6' cos'i/' + A sin-i/^) — 

h d — 
+p-^-^{Xcos6smil) + cosV')Pst- 



Xh' 



UiPi = p— sin 61 Pst, 



Eq. (15.61) becomes 



A rf2p , d 



2a de^ de 



X^ 
2a 



cot 9 — A (cos 9 + hes) sin 9 



P.t = 0. 



(5.9) 



(5.10) 



(5.11) 



The normalized solution of this equation assumes the form 



P^^(9) = Csin^exp[-2aiyefr 



(5.12) 



where Wcs{9) = Wcs{9) / mHa = (1/2) sin^ 9 — hcscos9 is the dimensionless effective energy 
of the nanoparticle and 

exp[a{l + hls)] 



C 



TT-^ erfi[v^(l + /icff)] + erfi[v^(l - h^g)] 



(5.13) 



is a normalizing constant derived from the condition 2it Jq d9Pst{9) = 1. Here erfi(2;) stands 
for the imaginary error function defined as erfi(2;) = {2/y/TT) /q^ (ixexp(x^). 

Finally, using the definition of the steady-state magnetization, poo = 27r Jq d9 cos9Pst{9), 
we obtain 

_ exp[a(l + /j-cff)^] - exp[a(l - h^sf] _ ~ 

V vraerfi[v/^(l + M]+erfi[v/^(l-M] 
Since erfi(2) = 2{z + /3 + ...)/ y/^ if z <^ 1 and erfi(2) = exp{z^){l/z + l/2z^ + . . .)/^ 

if 2; 3> 1, in the case of low anisotropy barrier (a ^ 1) this formula yields poo = 2a/iefr/3, 
and in the case of high anisotropy barrier (a ^ 1) and small effective field {\hes\ ^ 1) it is 
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reduced to the formula fioo = tanh(2a/iefr), which coincides with that derived from the MFPT 
approach. As an illustration of the accuracy of fl5.14p and the applicability of the MFPT 
approach for describing the magnetic relaxation in nanoparticle systems, we depict in Fig. 4 
the theoretical and the numerical results for the dependence of the induced magnetization 
/ioo on the anisotropy barrier height a = mHa/2kBT. 

C. Numerical verification 

As in the simulation of the MFPTs, described in the previous section, the Langevin equa- 
tions (12.101) were used to compute the mean magnetization at high frequencies. However, 
the numerical simulation in this case differs from the case studied before in two impor- 
tant aspects. Firstly, instead of an ensemble average over magnetic nanoparticles a single, 
stochastic magnetic moment was averaged over time, thereby making use of ergodicity. Sec- 
ondly, instead of using a reflecting and an absorbing boundary, here two reflecting boundaries 
were located aX 6 = O.OItt and at 6 = 0.997r. For fixed parameter values Q = 10, h = 1, 
p = — 1,A = 0.5 four trajectories were simulated for each value of a. Each trajectory was 
initialized with 6 = O.OSvr and ip = tt. We introduced two circles as marks on the sphere 
at 6* = 0.27r and 6 = O.Stt and defined a sign change of the magnetic moment as a crossing 
of the 6 = O.Stt [9 = 0.27r) mark, provided the magnetic moment had been in the 6 < 0.27r 
{6 > O.Svr) domain before. Each trajectory was run until two hundred such sign changes had 
occurred. The projection of the magnetic moment on the easy axis was summed over all time 
steps and divided by the number of time steps at the end of the simulation. The convergence 
of the mean magnetic moment along these trajectories was monitored and found to have 
converged after two hundred sign changes. From the four values for the mean magnetization 
a mean value and a standard deviation were determined. Fig. 4 depicts that the simulation 
and the analytic result of the steady-state theory at high frequencies, equation (I5.14p are 
in very good agreement, indicating that the high frequency limit is already obtained for 
n = 10. 
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VI. CONCLUSIONS 



We carried out a comprehensive study of the two-dimensional MFPT problem for the 
magnetic moment of a nanoparticle driven by a magnetic field rapidly rotating in the plane 
perpendicular to the easy axis of magnetization. Our approach is based on the equations 
(13. 7p and (13. 8p for the MFPTs that we derived from the backward Fokker-Planck equation in 
the rotating frame. In the high-frequency limit, we solved these equations analytically and 
calculated the MFPTs for the nanoparticle magnetic moment in the up and down states. 
The main finding is that a rapidly rotating field influences the MFPTs due to the change 
of the potential barrier between these states, which occurs under the action of the static 
effective magnetic field applied along the easy axis of magnetization. We showed that the 
magnetic field rotating in the clockwise (counter-clockwise) direction increases the MFPT 
for the magnetic moment in the up (down) state and decreases it for the magnetic moment 
in the down (up) state. Our theoretical predictions are in good agreement with the results 
obtained by numerical solution of the effective stochastic Landau-Lifshitz equations. 

In addition, we applied the derived MFPTs to study the features of magnetic relaxation in 
nanoparticle systems caused by a rotating magnetic field. We established that in the case of 
a large anisotropy barrier this field always decreases the relaxation time and magnetizes the 
nanoparticle system. The magnetization grows as the temperature decreases and its direction 
is uniquely determined by the direction of field rotation. Solving the forward Fokker-Planck 
equation in the case of a rapidly rotating field, we calculated also the magnetization for 
an arbitrary anisotropy barrier. The theoretical results are in excellent agreement with 
the numerical ones and confirm the applicability of the MFPT approach for describing the 
magnetic relaxation in systems of high-anisotropy nanoparticles driven by a rapidly rotating 
magnetic field. 

Due to the selective change of the noise-induced stability of a rapidly driven magnetic 
moment of a nanoparticle, as represented by corresponding mean first-passage times, the 
results herein can be used for potential applications in magnetic recording technology. The 
relative stability of the magnetic moment in the up and down states can be suitably con- 
trolled by either changing the temperature T of the environment (thereby changing the 
anisotropy barrier height) or upon varying the strength of a rapidly rotating magnetic field. 
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FIG. 1: Dependence of the ratio Tg-sep/'^'o- on the dimensionless anisotropy barrier height a for 
different values of the dimensionless amplitude h of the rotating magnetic field. The numerical 
calculations of the integrals in the relations (j4.7p and (|4.9p were carried out for $7 = 10, p = +1, 
= O.Ovr, = O.lvr, 9' = O.lvr if cr = +1, and 6' = O.Ovr if cr = —1. The broken curves (green 
and blue online) that cross the horizontal asymptote (black online) correspond to a = +1, and the 
broken curves (green and blue online) that lie below the asymptote correspond to a = —1. The 
solid curve (red online) represents the case h = for which T+iscp/^+i = T_isep/'?'-i- We note 
also that in all cases the ratio To-sep/'^'o- does not depend on A. 
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FIG. 2: The natural logarithm of the dimensionless MFPT T+i as a function of the parameter a. 
The thick curves represent the exact theoretical results obtained from (j4.7p . and the thin curves 
depict the approximate high barrier limit given by (j4.8p . The symbols indicate results from the 
numerical simulation of 4 x 10^ runs of the effective stochastic Landau-Lifshitz equations (|2.1U|) 
with the initial conditions ipo = 0, 9q = O.OSvr and with the absorbing boundary at = O.Svr. 
The broken curves and circular symbols (red online) correspond to = 5, and the solid curves 
and triangular symbols (blue online) correspond to 17 = 10. In accordance with the theoretical 
assumption, the analytical results approach to the numerical ones with increasing of the field 
frequency. 
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FIG. 3: The natural logarithm of ATg- as a function of the dimensionless amplitude h of the rotating 
magnetic field. The solid and broken curves represent the theoretical results obtained from the 
relation (j4.7p for a = +1 and a = —1, respectively. The values of the other parameters are a = 5, 
0, = 10, and p = +1. The symbols (in color online) depict the results obtained from the numerical 
simulation of Eqs. (j2.10p for different values of the damping parameter A. In full agreement with 
theoretical predictions, ln(AT(j) as a function of h decreases if o" = +1, increases if cr = —1, and 
does not depend on A. 
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FIG. 4: The dimensionless induced magnetization function of the anisotropy barrier height 

a for $7 = 10, h = I, p = —1, and A = 0.5. The sohd curve (red onhne) and the broken curve (blue 
onhne) represent the induced magnetization defined by the formulas (j5.14p and (15. 4p with (14. 7p . 
respectively. The triangular symbols indicate results obtained from the numerical simulation of 
Eqs. (j2.10p . As seen, the theoretical results that follow from the explicit solution of the Fokker- 
Planck equation are in excellent agreement with the numerical results. A small systematic shift 
of the induced magnetization derived within the rate theory arises from an overestimation of the 
absolute values of the average magnetic moment in the up and down states. 
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